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Abstract 

We investigated the equilibrium properties of a one-dimensional system of classical particles 
which interact in pairs through a bounded repulsive potential with a Gaussian shape. Notwith- 
standing the absence of a proper fluid-solid phase transition, we found that the system exhibits a 
complex behaviour, with "anomalies" in the density and in the thermodynamic response functions 
which closely recall those observed in bulk and confined liquid water. We also discuss the emer- 
gence in the cold fluid under compression of an unusual structural regime, characterized by density 
correlations reminiscent of the ordered arrangements found in clustered crystals. 
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I. INTRODUCTION 



Renewed interest has recently emerged in the thermodynamic and transport properties of 
one-dimensional (ID) systems, also in view of their potential relevance to the modelling of in- 
teracting colloidal particles diffusing within narrow channels [l| • The spontaneous emergence 
in such systems of (even intense) local ordering phenomena, not necessarily accompanied by 
proper phase transitions, is a topical issue that continues to be investigated also in relation 
to its effects on the dynamics of the particles 2j. In this regard, a landmark result still is the 
celebrated van Hove's theorem [3j according to which a system of identical particles with a 
hard core and pairwise interactions extending over a finite range has no phase transition in 
ID, in the absence of external fields. The general issue of the existence (or non-existence) of 
phase transitions in ID systems with short-ranged interactions has been critically reviewed 
in more recent years, among others, by Cuesta and Sanchez 4|. A way to potentially escape 
the negative verdict of van Hove's theorem is to allow for bounded repulsive interactions be- 
tween particles. The simplest of such models is a system of penetrable spheres (PS), whose 
interaction potential takes a finite positive value whenever two particles overlap, while van- 
ishing otherwise 5|, \&. In the absence of an infinite repulsive barrier, interactions are no 

n 

longer restricted to nearest neighbours. As a result, at variance with the Tonks gas [7|, no 
analytical solution of the PS model is currently known. Fantoni |8[ has recently provided 
robust arguments against the existence of a phase transition in a ID system of penetrable 
particles with a negative short-ranged square well outside the repulsive core jjj. 

A largely studied variant of bounded repulsive interactions is the Gaussian potential 10 ]. 
Its continuing interest also stems from its frequent applications as a model interaction in 
soft-matter systems, such as the dilute solutions of highly ramified polymers in a good 
solvent (star polymers), where the "effective" repulsion between the centres of mass of two 



polymer chains can be typically described through a Gaussian law The corresponding 
Gaussian-core model (GCM) can at most exhibit fluid-solid and, possibly, solid-solid phase 
transitions since a liquid- vapour transition is clearly excluded by the absence of an attractive 
component in the potential. In fact, at low enough temperature three-dimensional (3D) 
Gaussian particles crystallize, upon compression, in a cubic solid structure, with either a 
face-centred or a body-centred symmetry [12j. However, a further increase of the density 
does eventually lead to a reentrant melting of the solid. This phenomenon has been observed 
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also in two dimensions (2D) where, at variance with the 3D case, a narrow hexatic region 
appears between the isotropic fluid phase and the triangular solid phase all along the melting 



line [13 1 . 



As for the phase stability of the GCM in ID, no rigorous proof of the absence of a phase 
transition has been presented so far. However, Fantoni conjectured that his arguments 
against the existence of a phase transition in a ID system of attractive penetrable spheres 
apply equally well, under appropriate conditions, to a larger class of model fluids with a 
decaying long-ranged repulsive tail, including the GCM [8]. 

Whether or not the GCM fluid crystallizes in ID, one can safely expect a far from triv- 
ial phase behaviour because of the ever-increasing occurrence, under compression, of soft 
particle-core overlaps. A natural candidate for the macrostate of minimum Gibbs free en- 
ergy is a diffusely ordered arrangement formed by more or less equally spaced particles, as 
is the case of ID hard rods at high density [l4]. However, the thermodynamic competition 
between energy and entropy may also lead, in a fluid with a soft bounded potential, to even 
more complex arrangements such as those found in "clustered" or "tower" crystals, where 



two or more (superimposed) particles are confined within the same cell 11). 



II. MODEL AND METHOD 

We consider a system of point particles repelling each other, at a relative distance r, 
through a Gaussian pair potential 

u(r) = eexp [-{r/af] , (1) 

where e and a fix the energy and length scales, respectively. The GCM potential has an inflec- 
tion point at ro = <r/V2, where its curvature changes from concave to convex. Correspond- 
ingly, for r < r the strength of the force between two particles, f(r) = — dw(r)/dr, decreases 
as the particles approach each other. However, it is over a larger range (0 < r < cr) that the 
local virial function, rf(r), and thus the contribution of a pair of interacting molecules to 
the pressure of the system, decreases when the separation also decreases. This condition is 



typically taken as the signature of "core softening" [15|. As is well known, a core-softened 
- in the just specified sense - potential can generate a density anomaly, associated with 
a negative thermal expansion coefficient. This circumstance was originally verified for the 
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GCM fluid in 3D by Stillinger and Weber 16[. Indeed, a number of thermodynamic, struc- 
tural and dynamical properties of the GCM fluid have been already found to exhibit, both 
in 2D and 3D, "waterlike" anomalies which render the behaviour of this model qualitatively 
different from that of a "simple" (i.e., Lennard- Jones-like) fluid. 

To our knowledge, no systematic numerical study of the ID GCM fluid has been un- 
dertaken so far. In this paper we present the results of an investigation carried out with 
the Monte Carlo (MC) method in the isothermal-isobaric ensemble. The data we present 
were obtained with samples of N = 200 particles, unless otherwise specified; however, we 
also tested the stability and convergence of the results with larger samples of 500 and 1000 
particles. No qualitatively significant changes emerged upon increasing the "size" of the cal- 
culation. We collected data over trajectories of, typically, five million sweeps, every sweep 
consisting of N + 1 elementary MC moves, including one attempt, on average, to change 
the volume V of the system. During the equilibration runs, we adjusted the maximum par- 
ticle displacement and the volume change so as to maintain the acceptance rates of both 
types of moves close to 50% in the production runs. We carried out our simulations along 
a number of isothermal and isobaric paths, from low to high density and from high to low 
temperatures, continuing at each state point from the last system configuration produced 
in the previous run. We computed, for given values of T and P, the following properties: 
the average number density, n = N/ (V); the average energy, E = Nk B T/2 + (U), where k B 
is Boltzmann's constant and U is the total potential energy; the specific heat at constant 
pressure, Cp = T(ds/dT)p, where s is the entropy per particle; the isothermal compress- 
ibility, Kt = —v~ 1 (dv/dP)T, where v = n" 1 is the average volume per particle; the thermal 
expansion coefficient, ap = v~ l (dv/dT)p. The three thermodynamic response functions 
were obtained as thermal averages of covariances of the fluctuating variables (energy and 
volume): C P /k B = ([A(E + PV)} 2 ) / [k B T) 2 } a T /k B = (A(E + PV)AV) / [(V)(k B T)% and 
K T = ({AV) 2 )/ [(V)(k B T)}, where AX = X - (X). 

We also calculated the radial distribution function (RDF), g{\r\) = v(^2 k , 1 8(xk—Xi—r)), 
and the associated structure factor, S(\q\) = 1 + n dx exp(-iqx) [g(\x\) — 1]. 

The Monte Carlo study was complemented with the results obtained with the 
hypernetted-chain (HNC) approximation and with exact total-energy calculations carried 
out at T = with increasing pressure for several candidates for the solid phase, in order to 
gain insight into the preferred forms of particle aggregation at low temperatures. 
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In the following, we shall also make use of reduced units for density, temperature, and 
pressure: p = na, r = k^T/e, and II = Pa/e. 



III. RESULTS 

In the GCM fluid the growth of density correlations upon compression is eventually 
frustrated by the finite strength of the repulsion between particles and by the decreasing 
strength of their mutual force as they approach each other. We anticipate that for densities 
larger than the value corresponding to a nearest-neighbour (NN) distance roughly equal 
to 1.5cr (II ~ 0.35), a further increase of the pressure causes a suppression, rather than a 
sharpening, of the local structure of the fluid since in a sufficiently dense environment the 



overlapping of two particles entai 
penalty. Lang and coworkers 



s an entropy gain that is larger than the associated energy 
17] coined the term "infinite-density ideal gas" to represent 
the uncorrelated {g{r) = 1) behaviour of Gaussian particles at very high densities. 

Such a trend is already manifest in the RDF calculated through the HNC approximation 
at the reduced temperature r = 0.1 (see Fig. [1]) which, though quantitatively accurate only 
for p > 2, yields the correct qualitative trend of the local fluid structure with increasing 
pressure. We see from the picture that the NN distance, corresponding to the position 
of the first maximum in the RDF, steadily decreases upon compression but its statistical 
definition as a microscopic length scale is highest at an intermediate pressure corresponding 
to a reduced density p ~ 0.5. 

A. Phase stability at T = 

Before examining the MC findings, it is worth considering the thermodynamic behaviour 
of the fluid at zero temperature. Stillinger proved that in the limit of low T and n the 
thermodynamics of the GCM fluid reduces, for any space dimensionality, to that of hard 



spheres with a temperature-dependent diameter d(/3) = ay / hi(/3e), where /3 = (k-QT) -1 10]. 
In one dimension, the hard-rod behaviour sets in for pns — cr/d; below this threshold the 
equilibrium properties of the ID GCM fluid can be mapped onto those of the Tonks gas. 
However, pns is too small a density (about 0.40 for r = 0.002, corresponding to a pressure 
slightly smaller than 0.05) for the hard-sphere connection to be of any relevance to the 
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FIG. 1: Radial distribution function in the HNC approximation plotted as a function of the distance 
for a reduced temperature r = 0.1 and for reduced densities p = 0.1,0.2, . . . , 1.9,2 (the density 
increases along the direction indicated by the arrow). 

present inquiry into the existence of phase transitions in the ID GCM fluid at high pressure. 

The competition for thermodynamic stability between single-occupancy and clustered 
solids is ruled by the chemical potential which, at zero temperature, is equal to the enthalpy 
per particle: /i — e + Pv; moreover, any crystalline phase consists of just one microstate 
(modulo a translation). We first consider an arrangement of equally spaced particles: 

/i! = mini Vexp ( ~^-] + — \ ■ (2) 

We then compare fj,i with the chemical potentials of clustered phases with equidistant lumps 
of k — 2, 3, . . . particles, all lying at the same positions: 



k-1 



fc 2>p(-^H+- ■ ( 3 ) 

m=l x 7 ) 




FIG. 2: Excess chemical potential, relative to the value (/ii) in an ordered state of equally spaced 
particles, plotted as a function of pressure at T = 0, for clustered crystalline phases: continuous 
blue line, 2-clustered crystal; dashed cyanide line, 3-clustered crystal; long-dashed green line, 4- 
clustered crystal. The dotted red line yields the excess chemical potential of the state formed by 
regularly alternating singles and pairs. 

A further possibility that we have considered is a state in which isolated particles regularly 
alternate with pairs: 

1 . (4 ^ / 9m 2 \ 5 / 9m 2 \ p] 

/^ 2 ) = 3+mm - ^ exp^-— j+- ^ exp J + - . (4) 

L m=l,3,... v 7 m=2,4,... v 7 J 

In Fig. [2] we plot all such chemical potentials as a function of the pressure. The comparison 
shows that the crystal formed by equidistant particles is the most stable phase at zero 
temperature, whatever the pressure. No better solution was obtained upon allowing for a 
small fixed separation between the particles forming a pair in either the 2 — 2 — 2 — 2 — ... 
or the 1 — 2 — 1 — 2 — ... arrangement. This result is a strong indication of the absence 
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of any kind of phase transition in the ID GCM fluid. However, the free-energy penalties 
associated with the nucleation of clustered solids progressively decrease with the pressure, 
the faster the smaller the number of particles in a given cell. Hence, one cannot exclude a 
priori the possibility that, at non-zero temperature and high enough pressure, such phases 
may enter the thermodynamic game for entropic reasons. In passing, we note that for n > 3 
the 2-clustered crystal is almost degenerate with the single-occupancy crystal. 



B. Thermodynamic properties 

Following the computational procedure described in Sect. II, we investigated the phase 
behaviour of the ID GCM fluid at reduced temperatures 0.002 < r < 1 and for reduced 
pressures 0.05 < n < 3.5. We verified that neither the number density nor the energy exhibit 
any (even rounded-off) jump discontinuity over the explored domain. We also performed a 
chain of isothermal simulations at high pressure starting from an ordered 2-clustered con- 
figuration but did not find any sign of hysteresis in either the density or the energy. Hence, 
on a strictly thermodynamic basis, the equilibrium state of the investigated model is always 
that of a fluid. However, we found out that at low enough temperatures this fluid under- 
goes remarkable structural changes, with an anomalous (i.e., nonstandard) thermodynamic 
behaviour which, in many respects, closely resembles that observed in liquid water, when 
gradually cooled from above the freezing temperature into the deeply metastable region 
(possibly in a confined environment). 

The signature of such a "complex-fluid" behaviour is the emergence of a volumetric 
anomaly in the dense fluid: as shown in Fig. [3j at low pressure the density decreases mono- 
tonically with the temperature; however, for II ^ 1.2 the trend turns from decreasing to 
increasing over a limited temperature range T min (P) < T < T max (P). Correspondingly, a 
maximum shows up in n(T,P) at T = T max (P) which, for n < 1.5, is preceded by a min- 
imum at a nonzero temperature, T m i n (P). For n > 1.5, the density increases with T from 
T = 0, until it inverts its trend at the temperature of maximum density (TMD). 

A maximum inthe density of the GCM fluid had already been observed in higher space 
dimensions 



13 



16j | : in this respect, it is altogether remarkable to verify the persistence of 
this feature in one dimension. But even more remarkable is the emergence of a shallow 
minimum in n(T) over a range of pressures: the fluid behaves in an anomalous way for 
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FIG. 3: Reduced number density plotted as a function of the reduced temperature with increasing 
pressure: triangles (black), IT = 1; stars (blue), II = 1.2; squares (cyanide), II = 1.35; crosses 
(green), IT = 1.5; circles (magenta), II = 2; tripods (red), II = 2.5. To improve readability, the 
curves were shifted vertically by changing amounts; the separation between two major tick marks 
along the left vertical axis is equal to 0.01, while the absolute scale can be fixed through the value 
of each curve at r = 0.08, that is reported along the right vertical axis. 

T min (P) < T < T max (P) but turns "normal" again below the temperature of minimum 
density (TmD). In real liquids a density minimum has been detected in even fewer cases 
than those, already rare, where a maximum is observed. Two noteworthy examples are 



bulk tellurium [18|, also mixed with sulphur |l9|, and (supercooled) confined water [20h22| . 



where the emergence of a minimum has been associated with the formation of a defect-free 



"random tetrahedral network" 



23l |. However, the profound difference between tellurium 



and water (as noted by Angell, tellurium is not known as a network liquid 24|) makes it 



clear that even the TmD is caused by some more basic feature of the effective interaction 
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FIG. 4: Location of the density and response functions extrema in the (P,T) plane. Density: open 
red circles, temperature maxima (TMD); solid red circles, temperature minima (TmD). Coefficient 
of thermal expansion: magenta tripods, temperature minima; open magenta squares, pressure 
maxima; solid magenta squares, pressure minima. Isobaric specific heat: open blue triangles, 
temperature maxima; open blue squares, pressure maxima; solid blue triangles (inset), temperature 
minima. Isothermal compressibility: open black squares, temperature maxima; solid black squares, 
temperature minima. The locations of the symbols correspond to numerical estimates of the 
thermodynamic coordinates of the extrema obtained via Monte Carlo simulations along either 
isothermal or isobaric paths. The lines through the data are guides to the eye. 

potential. 

The relevant thermodynamic loci corresponding to the location of the extrema of the 
density and of the response functions are plotted in Fig. |H The volumetric anomaly region, 
corresponding to a negative value of the coefficient of thermal expansion, is the region 
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bounded by a maximum and (possibly) a minimum of the density; all along this boundary 
the expansivity vanishes. The TMD and TmD lines are seen to originate from point C (with 
coordinates lie — 1-169 and tq ~ 0.029), where the first and second isobaric temperature 
derivatives of the density both vanish [ctp = (dap/dT)p = p c = 0]. This point marks the 
onset of the volumetric anomaly in the fluid and is apparently reached with infinite slope 
along the TMD and TmD lines. The TMD line is seen to pass through a maximum (M), with 
coordinates 11m — 1.267 and tm — 0.037; correspondingly, the TMD locus has a positive 
slope for Pq < P < Pm, while decreasing asymptotically to zero with increasing pressure 
for P > Pm- Instead, the TmD shows a rapidly decreasing monotonic trend on compression 
and appears to vanish at a finite reduced pressure just larger than 1.5. 

The temperature Tm is the highest temperature at which the fluid manifests the volumet- 
ric anomaly, associated with the emergence of a single or double extremum in the density. 
This anomaly obviously reverberates in the thermodynamic response functions. Rolle's the- 
orem applied to an isothermal or isobaric "cut" of the anomalous region, producing two 
distinct intersections with the density extrema lines where ap = 0, implies that the expan- 
sivity must have a minimum at some point in between. The isobaric minimum of ap as 
a function of T (see Fig. [5]) goes along (though not exactly coinciding) with an inflection 
point of the density; both such features were found to survive well outside the anomalous 
region, down to II ~ 0.6 (see also Fig. [3]). The minimum-expansivity locus radiates out of 
this region at the point C of confluence of the TMD and TmD lines (see Fig. H]). 

Sastry and coworkers [25J] have showed that the isothermal compressibility increases 
on cooling at any point along the TMD line wherever its slope is negative and the 
second temperature derivative of the volume, evaluated at the same point, is positive 
[((9 2 u/9T 2 )p jat TMD = v(dap/dT)p^ at TMD > 0]. The above conditions apply for P > Pm- 
A similar thermodynamic constraint implies that the isothermal compressibility increases 
on heating at any point along the TmD line where the slope is negative and such also is 
the second temperature derivative of the volume, which is the case of the ID GCM fluid. 
The necessary consequence of this twofold constraint is the emergence of a maximum in 
the isothermal compressibility at some temperature in the interval T min (P) < T < T max (P). 
As shown in Figs. H] and El the maximum persists even for II > 1.5, when the minimum 
of the density has disappeared. We can see that, as originally predicted by Sastry and 
coworkers 25J, the locus of the compressibility maxima intersects the TMD line exactly at 
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FIG. 5: Thermal expansion coefficient plotted as a function of the reduced temperature for different 
pressures (left panel) and of the pressure for different temperatures (right panel); left panel: same 
legends as in Fig. [3l right panel: triangles (black), r = 0.002; stars (blue), r = 0.005; squares 
(cianide), r = 0.03; crosses (green), r = 0.05; circles (magenta), r = 0.1; tripods (red), r = 1. 

its maximum M, where (dT max (P)/dP)p== 0. 



Thermodynamic consistency implies 



25j that the anomalous temperature dependence of 



Kt should also reflect in the behaviour of Cp as a function of T at fixed P (see the left panel 
of Fig. EJ). As the temperature drops to zero, the isochoric specific heat approaches unity 
(i.e., the value of an interacting ID fluid at T = 0), as also does the isobaric specific heat 
since Cp = Cy + Tvap/K T (see Fig. [5}; Cp shows then a maximum, whose (P,T) locus 
runs very close to the isobaric minimum-expansivity line (see Fig. H]). We also note that the 
height of the Cp maximum changes very little with either the pressure or the temperature, 
an aspect that we shall return to soon. 
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FIG. 6: Isothermal compressibility plotted as a function of the reduced temperature for different 
pressures (left panel) and of the pressure for different temperatures (right panel); same legends as 
in Fig. [5j To improve readability, the curves in the left panel were shifted vertically by changing 
amounts; the separation between two major tick marks is 0.05, while the absolute scale can be 
fixed through the value of each curve at r = 0.1, that is reported along the right vertical axis. 

Upon heating the system, the thermodynamic properties of the fluid should asymptot- 
ically approach their ideal-gas values. As a result, in order for their low-temperature be- 
haviour to match the high-temperature regime, a second extremum must develop in all of 
the three response functions, besides the one at low temperature that we have already com- 
mented on: this second extremum is a maximum in otp which, as T — > oo, must tend to zero 
as 1/T; instead, a minimum is found in both Kt and Cp which should approach 1/P and 
3/2, respectively. Hence, the second extremum exhibited by the response functions is the 
necessary outcome of the existence of two different limiting behaviours and, as such, does 
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FIG. 7: Isobaric specific heat plotted as a function of the reduced temperature for different pressures 
(left panel) and of the pressure for different temperatures (right panel) ; same legends as in Fig. [3 

not arise from any additional independent anomaly of the fluid at high temperature. 

The line of the isobaric temperature minima of Kt has also been traced in Fig. HI where 
it is apparent that this locus and the locus of the isobaric temperature maxima of Kt stem 
out of the same point, located at Hr t — 0.84 and t Kt ~ 0.11. These two lines get closer 
again at high pressure. 

The right panels of Figs. El El El and [S] show how the fluid reacts to an isothermal com- 
pression. Both the density and the energy were found to steadily increase with increasing 
pressure. Moreover, the isothermal compressibility is seen to decrease monotonically with P, 
its value being not significantly affected by the temperature over the explored range. Instead, 
ap and Cp exhibit a more complex behaviour. As for the expansivity, the thermodynamic 
relation (dap/dP)p = —(dKT/dT)p implies that isobaric temperature extrema of the com- 
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FIG. 8: Isochoric specific heat plotted as a function of the reduced temperature for different 
pressures (left panel) and of the pressure for different temperatures (right panel); same legends as 
in Fig. [5j 



pressibility should map onto isothermal pressure extrema of the expansivity [26[. In fact, 
Fig. H] shows that the isothermal pressure minima (maxima) of ctp do actually fall, to within 
the numerical uncertainty of the calculations, along the locus of the isobaric temperature 
maxima (minima) of Kt- 

A maximum is also present, at low temperature, in the isobaric specific heat as a function 
of P at constant T (see the right panel of Fig. [7]); as expected, this quantity is 3/2 at P = 
and, correspondingly, Cy = 1/2. We found that the locus of the pressure extrema of Cp runs 
on top of the temperature minima locus of the expansivity (see Fig. H]). The thermodynamic 
relation 



dC P 
dP 



-vT 



ap + 



dap 



(5) 
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states that, if ap = 0, an isobaric temperature extremum of ap should coincide with an 
isothermal pressure extremum of Cp. This is precisely what happens at the confluence point 
of the TMD and TmD lines. Elsewhere, the balance between ap and a negative temperature 
derivative of ap locates the zero of (dCp/dP)T at lower temperatures with respect to the 
locus of vanishing (dap/dT)p. However, upon moving inside the anomaly region, the first 
term on the r.h.s of Eq. |5]turns out to be very small (of the order of 10~ 4 ) for pressures just 
greater than lie, and becomes even smaller with increasing pressure because of the gradual 
drop of both the volume and temperature. As a result, the pressure maxima locus of Cp 
runs very close to the temperature extrema loci of ap and of Cp itself (see Fig. H]). This 
latter circumstance, i.e., the nearby coincidence (to within the uncertainty of the present 
calculations) of the pressure and temperature maxima loci of Cp, can be explained by the 
implicit function theorem: 



where X(P,T) is a generic property of the system and the second partial derivative on the 
r.h.s. of Eq.(j6]) is evaluated along a constant-X thermodynamic path. In general, whenever 
X has an extremum as a function of P at given T (i.e., (dX/dP) T = 0), the derivative 
(dT/dP) x , evaluated along a constant-X path, vanishes as well (a similar argument natu- 
rally holds upon inverting P with T). This obviously implies that a pressure (temperature) 
extremum of X is not necessarily associated with a temperature (pressure) extremum of 
the same quantity. However, it can be shown that, when X(T, P) displays a "ridge" of 
equal-height extrema, both the pressure and the temperature derivatives of X vanish there 
(and vice versa), which means that the extremum is such along both thermodynamic axes. 
To a very good approximation, this is actually what happens for the isobaric specific heat 
of the ID GCM fluid. 

At very low temperature the peak of the isobaric specific heat is preceded by a broad 
minimum, actually almost a plateau at one. This latter circumstance is explained by the 
very small values attained by the quantity vTap/ K?, which yields the difference between 
Cp and Cy; moreover, ap vanishes at a point located within the implicated pressure interval. 
In passing, we note that Cp = Cy along the TMD and TmD lines. 

One more aspect of this thermodynamic scenario still needs to be outlined. The volu- 
metric anomaly is accompanied by an anomaly of the entropy in that, whenever ap < 0, the 




(6) 
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entropy increases with P at constant temperature since (dS/dP)r = —(dV/dT)p. In order 
to illustrate this aspect, we calculated the total entropy per particle of the fluid through the 
Euler relation, s = j3(u + Pv — fj). As for the chemical potential, given its value at some 
reference state, one can calculate its value at any other state upon integrating a suitable 
thermodynamic property along either an isothermal or isobaric path, under the condition 
that no coexistence locus is being crossed along the path: 

MT, P 2 ) = fj,(T, Pi) - / 2 dP v{T, P) , (7) 



fi(T 2 ,P) _ t i(T 1 ,P) 
T 2 Ti 



Pi 



dT 



Ti 



u(T } P) + Pv(T,P) 

J^2 



(8) 



Equations [7| and [8] readily follow from the Euler and Gibbs-Duhem relations. The chemical 
potential at the reference thermodynamic state (a dilute-fluid state) can be estimated using 



Widom's particle-insertion method |28[. Upon performing a NVT simulation for p = 0.2 
and r = 0.1, we obtained the value = —0.7341 at a reduced pressure of 0.0296. Using 
this value as a parameter in a NPT simulation, carried out at the same temperature as in 
the const ant- volume simulation, we obtained the value (3fi = —0.7348, which coincides with 
the corresponding NVT estimate to within the statistical error of the calculation. 

The left panel of Fig. [9] shows the results for the total entropy per particle of the ID 
GCM fluid. As expected, at high temperature this quantity decreases with P, because 
of the increasing strength of the positional correlations, a process which typically leads to 
a reduction of the number of microstates available to the system in a given macrostate. 
However, as soon as the fluid enters the anomalous region, viz., when the temperature drops 
below Tm, the total entropy develops a minimum that is followed by a maximum at larger 
pressure, a state beyond which the entropy starts decreasing again. This counterintuitive 
behaviour is a direct consequence of the volumetric anomaly; correspondingly, the locations 
of the two pressure extrema change with the temperature, following the TMD and TmD 
lines. 

The right panel of Fig. [9] shows that a minimum is also present, at lower pressure, in the 
excess entropy of the fluid, S cx = S — S^, where is the entropy of an ideal gas having 
the same density and temperature of the GCM fluid. However, three notable differences 
show up in the thermodynamic behaviours of the excess and total entropy, respectively: i) 

(s ) 

the minimum of the excess entropy is located at Tl m tn — 0-35 and its position does not 
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FIG. 9: Total entropy (left panel) and excess entropy (right panel) per particle plotted as a function 
of the pressure for different reduced temperatures: triangles (black), r = 0.002; stars (blue), 
r = 0.005; squares (cianide), r = 0.03; crosses (green), r = 0.05; circles (magenta), r = 0.1. The 

1 /2 

total entropy reported on the left panel does not include the additive constant In [27rm/(e/i 2 )] ; 
the arrows indicate the positions of the extrema; a dotted line has been traced for II = 0.35 in the 
right panel. 

appreciably shift with T over the explored range; ii) for II > IT^ the excess entropy 
rises to zero and does not exhibit any other extremum; iii) the minimum persists in s ex 
even for T > T M . All such circumstances suggest that the pressure threshold at ^ 
may be a significant indication of the existence of two markedly different structural regimes 
in the fluid. We shall come back to this point in the following section. However, before 
concluding this outline of the thermodynamic properties of the model, we note that the 
excess entropy of the GCM fluid displays a minimum also in upper dimensions (2D, 3D). 
For temperatures slightly above the maximum melting temperature, this minimum occurs 
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at a reduced density corresponding to an average interparticle distance of approximately 
1.53a, a value that is fairly congruent with the values in ID, which range between 1.46a 
and 1.56o" as the reduced temperature rises from 0.002 to 0.1. In 2D and 3D the excess- 
entropy minimum falls just beyond the density corresponding to the maximum melting 
temperature, which unambiguously marks the border between two different thermodynamic 
conditions: a low-density regime in which the system behaves as a normal fluid and, upon 
compression, eventually freezes; a high-pressure regime in which, instead, the crystal re- 
melts, with increasing pressure, into a fluid whose properties are very different from the 
ordinary ones. 

C. Structural properties 

In this section we shall focus on the structural properties of the ID GCM at low tempera- 
ture (r = 0.002). Upon isothermally compressing the initially dilute fluid, the local order is 
more and more enhanced, and also extends to larger and larger distances. As seen in Fig. [TUJ 
the RDF already looks highly structured at a relatively low pressure (II = 0.05, p = 0.47), 
showing a quasi-crystalline profile characterized by a series of sharply defined coordination 
shells (each subtending an integrated conditional number density of one particle), with pe- 
riod equal to the average interparticle distance and no appreciable overlap between adjacent 
shells up to fairly large distances. Such a regularly equispaced arrangement is analogous to 
that exhibited by a ID gas of hard particles for packing fractions approximately larger than 
83% [14], a value corresponding to a reduced density of the present model equal to 0.33 at a 
reduced temperature of 0.002. In this respect, the core-softened GCM fluid matches and ex- 
tends to higher densities the quasi-crystalline phase behaviour observed in the "isomorphic" 
hard-core Tonks gas approaching close packing. 

As expected for an ordinary simple fluid, the peaks of the RDF get, with increasing 
pressure, sharper and taller; contextually, the size (Ro) of the domain over which the minima 
between adjacent peaks have, to practical effects, dropped to zero also expands with P. Such 
a trend comes to a stop and is eventually reversed when the pressure increases to values 
close to the minimum threshold discussed above for the excess entropy (see Figs. ITTland [T2|) . 
In fact, for II ~ 0.30, both the height of the first RDF peak as well as Rq attain a maximum 
(of about 11 and 20t>, respectively) and start decreasing thereon with P. Figure [TBI shows 
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FIG. 10: Radial distribution function at r = 0.002 plotted as a function of the distance between 
particles, relative to their average separation at the corresponding number density, for II = 0.05 
(blue) and II = 0.35 (red); the lower panel shows the decay of the total correlation function at 
larger distances on a semilogarithmic scale. 

that a similar inversion of the low-density trend is also observed in the other maxima of the 
RDF at reduced pressures close to 0.35. Hence, this threshold corresponds to the maximum 
growth and expansion of the local (hard-core-like) crystalline order in the fluid: in fact, for 
larger pressures, we observe a global attenuation of the density correlations which reduces 
the entropic distance of the fluid from its ideal-gas counterpart. 

The modifications of the local density profile outlined above prelude the emergence, pro- 
moted by the pressure, of a different spatial organization of the fluid. Figure [13] shows that 
the larger the distance of a given coordination shell from the particle sitting at the origin, 
the slower the decay of the corresponding peak height with increasing pressure. As a result, 
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FIG. 11: Radial distribution function at r = 0.002 plotted as a function of the distance between 
particles, relative to their average separation at the corresponding number density, for II = 1.5 
(blue) and II = 2 (red); the lower panel shows the decay of the total correlation function at larger 
distances on a semilogarithmic scale. 



the second maximum in g(r) eventually overcomes the first peak for reduced pressures ap- 
proximately larger than 1.4, i.e., as soon as the fluid enters the volumetric anomaly region. 
An analogous crossover between the fourth and third peak takes place for IT ~ 1.75. A 
novel type of local order gradually builds up in the fluid, which, for IT > 2.5, shows up in 
the RDF as the superposition of two distinct modulations, sustained by the odd- and even- 
numbered peaks, with spatial periods both equal to 2v but with different envelopes. Overall, 
the structure of the fluid is dampened by the pressure, but not everywhere in the same way: 
odd-numbered peaks get progressively lower and broader than even peaks; moreover, the 
heights of the even peaks are seen to decay monotonically with the distance, at variance with 
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FIG. 12: Radial distribution function at r = 0.002 plotted as a function of the distance between 
particles, relative to their average separation at the corresponding number density, for II = 2.7 
(blue) and II = 3.5 (red); the lower panel shows the decay of the total correlation function at larger 
distances on a semilogarithmic scale. 

those of the odd-numbered peaks which display a maximum at some intermediate distance 
R c ; this distance grows with the pressure and appears to saturate, for IT > 3, to a size 
corresponding to the first 12 coordination shells. For r > R c , the two trains of oscillations 
do eventually merge into one single oscillation with period v. 

An indication on the nature of the new spatial organization spawned by the fluid upon 
compression emerges from a typical snapshot taken at the reduced pressure II = 3: the 
presence of patches of particle pairs ("dimers"), regularly spaced with period 2v and width 
of the order of R c , is rather evident in Fig. [TH The fact that the even-numbered peaks of the 
RDF are taller and sharper indicates that the average distance between successive dimers 
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FIG. 13: Heights of the first five maxima of the RDF, identified with Roman numerals, plotted as 
a function of the reduced pressure at r = 0.002. 

is rather well defined as is the next-nearest-neighbour distance between "isolated" particles. 
On the other hand, the nearest-neighbour coordination shell looks less resolved because, if 
the particle located at the origin belongs to a dimer, one of its two first neighbours will 
be closer on average than the other one, belonging to an adjacent dimer. The onset of 
the pairing phenomenon at high pressure is further witnessed by the value of the RDF at 
r = 0, which gives a measure of the fraction of microstates with superimposed particles: 
this quantity vanishes (to more than six significant figures) for reduced pressures lower than 
2.3 but starts increasing monotonically thereon, being approximately 0.1 for II = 3. The 
crossover from a single-occupancy fluid to a dimer-rich fluid is continuous rather than sharp, 
i.e., it does not entail any thermodynamic transition. We obviously expect that "trimers" 
also form at still higher densities. 
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FIG. 14: A snapshot of part of the system at r = 0.002 and II = 3 (p = 1.8421); particles are 
represented as black dots and distances are relative to the average separation at the given density: 
we can clearly observe the prevailing presence of regularly spaced patches of "dimers" interspersed 
with ordered sequences of single particles. 

We conclude our discussion of the changes observed in the RDF upon compression with 
a note on the decay of this function at large distances, which provides an indication on 
the spatial extent of coherent density fluctuations in the fluid. To this end, we successfully 
fitted the envelope of more distant maxima to an exponentially decaying function. This 
very circumstance further confirms the fluid nature of the system even at low temperature 
and high pressure. The fit allowed us to extract the positional correlation length £, whose 
behaviour is shown as a function of P in Fig. [15j When plotted in units of the average 
separation between particles at a given density, this quantity exhibits a rapid growth for 
P < 0.35, the structural threshold where the tendency towards a single-occupancy solid- 
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FIG. 15: Correlation length, relative to the average separation between particles, plotted as a 
function of the reduced pressure at r = 0.002: dotted blue curve, 200 particles; continuous red 
curve, 500 particles. 

like arrangement produces its strongest effects. For larger pressures, after a slight decrease, 
the correlation length starts growing again, more slowly than at low pressures, with an 
approximately linear growth rate for the larger size investigated. We also verified that £ never 
exceeds half of the size of the simulated sample, up to the highest pressures investigated. 

The structure factor provides a complementary view of the density correlations in the 
fluid (see Fig. [16]) . At low pressure, the main peaks of S(q) are found at values of the scaled 
wavenumber qv ~ ± 2tt, corresponding to a large-distance oscillation of the RDF with period 
v. However, with increasing pressure new satellite peaks - centred at odd multiples of 7r 
- sequentially emerge on the high- wavenumber side, the symmetric peaks at ± 7r being the 
last to appear (for IT > 1.4, again approximately when the fluid enters the anomaly region). 
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FIG. 16: Structure factor plotted as a function of the reduced wavevector at r = 0.002 for II = 1 
(black), 1.5 (blue), 2 (cyanide), 2.5 (green), 3 (magenta), 3.5 (red); a specular spectrum is found 
for negative values of the wavevector. The dotted vertical lines lie in correspondence to the first 
four multiples of n. 

Such additional maxima initially concur to modify the profile of the first coordination shell 
in the RDF, the more so the larger the pressure. However, the increasing relevance of the 
peaks at ± tt, which eventually overcome in height the peaks at ± 3n for reduced pressures 
close to 3, also signals the emergence of an additional spatial modulation in the local density 
profile, with period 2v. 
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IV. CONCLUDING REMARKS 



In this paper we have illustrated the equilibrium properties of a model system of particles 
interacting, in one dimension, through a bounded repulsive potential with a Gaussian shape. 
The study of the model was largely carried out using numerical sampling techniques based 
on the Monte Carlo method. 

As far as the thermodynamic phase stability of the model is concerned, no indication 
emerged from our analysis of any singular behaviour possibly associated with a fluid-solid or 
solid-solid transition. However, this very circumstance allowed us to carry out a thorough 
study of the unexpectedly complex behaviour of the Gaussian-core model fluid in one di- 
mension, all the way from high to low temperatures, without incurring any "break" caused 
by phase changes. As a result, we could highlight the full thermal unfolding of a volumet- 
ric anomaly - i.e., the increase of the density upon isobarically heating the system - as a 
function of pressure, with the emergence of both a maximum and a minimum of the density 
with decreasing temperature. While the existence of a maximum in the density of the model 
had already been observed both in three and two dimensions, evidence of the existence of 
a minimum, at lower temperature and under appropriate pressure conditions, has never 
been reported and is being provided here, for the first time, for the one- dimensional fluid. 
This feature is more unusual than the maximum itself and has been recently observed in 
metastable states of real substances, most notably in confined supercooled water. 

In the light of such findings, we can safely conclude that a soft bounded repulsion like that 
modelled with a bare Gaussian potential, devoid of any additional repulsive or attractive 
component, possesses the "minimal" requisites that make it possible for the fluid to exhibit 
waterlike anomalies even in a one-dimensional hosting space. Sadr-Lahijany and cowork- 
ers 29(| had already described the nonstandard thermodynamic properties that emerge from 



a core-softened potential, explicitly designed to mimic the effect of hydrogen bonding: a re- 
pulsion with two length scales set by a hard core and by a finite negative shoulder, followed 
at larger distances by an attractive well. In this respect, the main result of this paper is that 
definitely much less is needed to "turn on" a volumetric anomaly, even in a one-dimensional 
fluid: in fact, a finite softened repulsion, modelled as a one-scale potential with a downward 
concavity at short distances and an upward one at larger distances, is sufficient to generate 
the density anomaly, also bounded by two temperature extrema, as well as a cascade of 
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related peculiarities in all of the three thermodynamic response functions. 

Crucial for the thermodynamic onset of the volumetric anomaly in the currently inves- 
tigated model is the average separation between particles at a given pressure: in fact, we 
may conjecture that, if such a distance is smaller than the distance r = cx/y/2 where the 
concavity of the Gaussian potential changes from upward to downward, the density will grow 
upon isobarically heating the fluid even at zero temperature, because the particles will find 
it (thermodynamically) more advantageous to reduce their average mutual distance, since 
this will ultimately entail a significant gain in the entropy of the fluid which overcomes the 
corresponding increase of internal energy. Such a "marginal" condition corresponds to a 
threshold density po — which, in the ordered ground state formed by equally spaced 
particles, is attained with a reduced pressure Ilo — 1.77. Indeed, a similar mechanism is at 
work in liquid water, in which the breaking of hydrogen bonds, with their hybrid attractive 
and repulsive nature [30j , produces a gradual "collapse" of nearest-neighbour molecules in- 
side the inner region of close-contact distances, with a corresponding increase of the density 
upon heating. 

On the other hand, we can also presume that, if the average separation between particles 
is larger than r , the density of the fluid will initially decrease upon heating the system from 
a state at T = (because the energy decreases as well) until the particles acquire enough 
kinetic energy to sample the inner region of the potential and thus discover that they may 
exploit there a more favourable condition, leading to an inversion of the density trend as 
a function of T. On this basis, we may then expect that for II < Ilo, and for not too low 
pressures, the anomalous region is bounded from below by a nonzero temperature at which 
the density shows a minimum. 

The investigation of the structural properties of the model also revealed an interesting and 
unusual scenario. Two structural regimes were distinguished in the equilibrium behaviour of 
the ID GCM: for reduced pressures lower than 0.35, the system behaves as a "normal" fluid 
in that the local order is progressively enhanced by an isothermal compression. Instead, 
upon trespassing the above threshold, a further increase of the pressure induces an overall 
attenuation of density correlations, an effect associated with the bounded nature of the 
repulsion. However, the ensuing approach of the system to the thermodynamic condition 
of an "infinite-density ideal gas" is accompanied by the emergence of a new type of local 
order, with pairs of almost superimposed particles giving rise to extended quasi- crystalline 
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clusters of "dimers" . This phenomenon is clearly resolved in the radial distribution function 
as an extra modulation, which is also signalled by the appearance of satellite peaks in the 
structure factor at relative wavenumbers equal to odd multiples of n. 
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